Volcano plot
diff_percent_isoform_volcano_padj <-
diff_percent_isoform_pval |>
volcano_plot(
x = diff_percent_isoform, y = p_adj, label = padj_label, FC_threshold = 1,
xlab = '∆ (% isoform)\n (KD - control)',
ylab = paste('-log10 (p.adj)')
)
diff_percent_isoform_volcano_padj |>
ggsave_multiple_formats(
outdir = figdir, height = 9, width = 9, fontsize = 7
)
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
diff_percent_isoform_volcano_p <-
diff_percent_isoform_pval |>
volcano_plot(
x = diff_percent_isoform, y = p, label = p_label, FC_threshold = 1,
xlab = '∆ (% isoform)\n (KD - control)'
)
diff_percent_isoform_volcano_p |>
ggsave_multiple_formats(
outdir = figdir, height = 18, width = 18, fontsize = 7
)
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27216 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1430 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27216 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: ggrepel: 1431 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27216 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1430 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27216 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: ggrepel: 1431 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27216 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1435 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
log2FC_percent_isoform_volcano_padj <-
diff_percent_isoform_pval |>
volcano_plot(
x = log2FC_percent_isoform, y = p_adj, label = padj_label,
repel_max_overlaps = 20,
FC_threshold = 1,
xlim = c(-6,6), x_num_breaks = 7
)
log2FC_percent_isoform_volcano_padj |>
ggsave_multiple_formats(
outdir = figdir, width = 9, height = 9, fontsize = 7
)
## Warning: Removed 2132 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2132 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2132 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2132 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 2132 rows containing missing values (`geom_point()`).
## Warning: Removed 28659 rows containing missing values (`geom_text_repel()`).
log2FC_percent_isoform_volcano_p <-
diff_percent_isoform_pval |>
volcano_plot(
x = log2FC_percent_isoform, y = p, label = p_label, FC_threshold = 1.2, #xlim = c(-5, 5),
repel_max_overlaps = 20
)
log2FC_percent_isoform_volcano_p |>
ggsave_multiple_formats(
outdir = figdir, width = 9, height = 9, fontsize = 7
)
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27523 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1135 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27523 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1135 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27523 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1135 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27523 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1135 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 2128 rows containing missing values (`geom_point()`).
## Warning: Removed 27523 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 1121 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
diff_percent_isoform_m3Cgenes_volcano_padj <-
diff_percent_isoform_m3Cgenes_pval |>
volcano_plot(
x = diff_percent_isoform, y = p_adj, label = padj_label, FC_threshold = 1,
xlab = '∆ (% isoform)\n (KD - control)',
ylab = paste('-log10 (p.adj)')
)
diff_percent_isoform_m3Cgenes_volcano_padj |>
ggsave_multiple_formats(
outdir = figdir, height = 9, width = 9, fontsize = 7
)
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
diff_percent_isoform_m3Cgenes_volcano_p <-
diff_percent_isoform_m3Cgenes_pval |>
volcano_plot(
x = diff_percent_isoform, y = p, label = p_label, FC_threshold = 1, xlim = c(-20, 20),
xlab = '∆ (% isoform)\n (KD - control)', repel_max_overlaps = 20
)
diff_percent_isoform_m3Cgenes_volcano_p |>
ggsave_multiple_formats(
outdir = figdir, height = 9, width = 9, fontsize = 7
)
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 370 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 48 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 370 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: ggrepel: 48 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 370 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 48 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 370 rows containing missing values (`geom_text_repel()`).
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: ggrepel: 48 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <e2>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <88>
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## conversion failure on '∆ (% isoform)' in 'mbcsToSbcs': dot substituted for <86>
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 370 rows containing missing values (`geom_text_repel()`).
## Warning: ggrepel: 43 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
log2FC_percent_isoform_m3Cgenes_volcano_p <-
diff_percent_isoform_m3Cgenes_pval |>
volcano_plot(
x = log2FC_percent_isoform, y = p, label = p_label,
FC_threshold = 1.5, p_threshold = .01,
xlim = c(-6, 6), x_num_breaks = 7,
repel_max_overlaps = 20
)
log2FC_percent_isoform_m3Cgenes_volcano_p |>
ggsave_multiple_formats(
outdir = figdir, width = 9, height = 9, fontsize = 7
)
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
## Warning: Removed 14 rows containing missing values (`geom_point()`).
## Warning: Removed 420 rows containing missing values (`geom_text_repel()`).
barplot
barplot_percentage_isoforms('S100A4')
## # A tibble: 54 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 Cont D N1
## 2 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 Cont D N2
## 3 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 Cont D N3
## 4 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… G N1
## 5 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… G N2
## 6 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… G N3
## 7 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… I N1
## 8 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… I N2
## 9 ENST00000354332.8 S100A4-201 ENSG0000019615… S100A4 siME… I N3
## 10 ENST00000368715.5 S100A4-203 ENSG0000019615… S100A4 Cont D N1
## # ℹ 44 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 12 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000368… perc… D G 3 3 -21.3 3.12 1.73e-4 3.46e-4
## 2 ENST00000368… perc… D I 3 3 -23.0 2.08 2 e-3 2 e-3
## 3 ENST00000481… perc… D G 3 3 -1.98 4.00 1.19e-1 1.19e-1
## 4 ENST00000481… perc… D I 3 3 -5.57 3.45 8 e-3 1.6 e-2
## 5 ENST00000354… perc… D G 3 3 42.8 3.04 2.51e-5 2.51e-5
## 6 ENST00000354… perc… D I 3 3 36.8 4.00 3.28e-6 6.56e-6
## 7 ENST00000468… perc… D G 3 3 6.96 3.72 3 e-3 3 e-3
## 8 ENST00000468… perc… D I 3 3 -16.5 2.63 1 e-3 2 e-3
## 9 ENST00000368… perc… D G 3 3 -3.91 2.35 4.6 e-2 9.1 e-2
## 10 ENST00000368… perc… D I 3 3 -0.199 2.15 8.59e-1 8.59e-1
## 11 ESPRESSO:chr… perc… D G 3 3 12.6 2.27 4 e-3 4 e-3
## 12 ESPRESSO:chr… perc… D I 3 3 15.5 2.47 2 e-3 3 e-3
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('RPS24')
## # A tibble: 18 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 Cont D N1
## 2 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 Cont D N2
## 3 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 Cont D N3
## 4 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… G N1
## 5 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… G N2
## 6 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… G N3
## 7 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… I N1
## 8 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… I N2
## 9 ENST00000372360.9 RPS24-202 ENSG0000013832… RPS24 siME… I N3
## 10 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 Cont D N1
## 11 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 Cont D N2
## 12 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 Cont D N3
## 13 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… G N1
## 14 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… G N2
## 15 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… G N3
## 16 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… I N1
## 17 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… I N2
## 18 ENST00000613865.5 RPS24-214 ENSG0000013832… RPS24 siME… I N3
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST000006138… perc… D G 3 3 26.7 3.99 1.20e-5 2.40e-5
## 2 ENST000006138… perc… D I 3 3 5.68 2.92 1.2 e-2 1.2 e-2
## 3 ENST000003723… perc… D G 3 3 -32.6 2.41 3.02e-4 6.04e-4
## 4 ENST000003723… perc… D I 3 3 -9.95 2.27 6 e-3 6 e-3
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('RPS14')
## # A tibble: 27 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 Cont D N1
## 2 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 Cont D N2
## 3 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 Cont D N3
## 4 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… G N1
## 5 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… G N2
## 6 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… G N3
## 7 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… I N1
## 8 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… I N2
## 9 ENST00000312037.6 RPS14-201 ENSG0000016458… RPS14 siME… I N3
## 10 ENST00000407193.7 RPS14-203 ENSG0000016458… RPS14 Cont D N1
## # ℹ 17 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000407193.7 perce… D G 3 3 -5.22 3.73 0.008 0.008
## 2 ENST00000407193.7 perce… D I 3 3 -7.12 3.94 0.002 0.004
## 3 ENST00000519690.1 perce… D G 3 3 0.353 3.21 0.746 1
## 4 ENST00000519690.1 perce… D I 3 3 -0.717 4.00 0.513 1
## 5 ENST00000312037.6 perce… D G 3 3 3.63 3.10 0.034 0.045
## 6 ENST00000312037.6 perce… D I 3 3 5.12 2.48 0.023 0.045
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('MDK')
## # A tibble: 45 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000395565.5 MDK-202 ENSG0000011049… MDK Cont D N1
## 2 ENST00000395565.5 MDK-202 ENSG0000011049… MDK Cont D N2
## 3 ENST00000395565.5 MDK-202 ENSG0000011049… MDK Cont D N3
## 4 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… G N1
## 5 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… G N2
## 6 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… G N3
## 7 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… I N1
## 8 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… I N2
## 9 ENST00000395565.5 MDK-202 ENSG0000011049… MDK siME… I N3
## 10 ENST00000395566.9 MDK-203 ENSG0000011049… MDK Cont D N1
## # ℹ 35 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000395566.9 perc… D G 3 3 -3.16 3.95 0.035 0.07
## 2 ENST00000395566.9 perc… D I 3 3 3.80 2.13 0.056 0.07
## 3 ENST00000407067.1 perc… D G 3 3 1.49 2.71 0.242 0.242
## 4 ENST00000407067.1 perc… D I 3 3 -7.63 3.47 0.003 0.005
## 5 ENST00000395565.5 perc… D G 3 3 3.88 3.80 0.02 0.039
## 6 ENST00000395565.5 perc… D I 3 3 0.116 3.81 0.913 0.913
## 7 ENST00000405308.6 perc… D G 3 3 2.93 3.82 0.045 0.063
## 8 ENST00000405308.6 perc… D I 3 3 3.29 3.90 0.032 0.063
## 9 ESPRESSO:chr11:3… perc… D G 3 3 -0.393 3.90 0.715 1
## 10 ESPRESSO:chr11:3… perc… D I 3 3 -0.357 2.22 0.752 1
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('SH3BGRL3')
## # A tibble: 27 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 Cont D N1
## 2 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 Cont D N2
## 3 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 Cont D N3
## 4 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… G N1
## 5 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… G N2
## 6 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… G N3
## 7 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… I N1
## 8 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… I N2
## 9 ENST00000270792.10 SH3BGRL3-201 ENSG000001426… SH3BGRL3 siME… I N3
## 10 ENST00000319041.6 SH3BGRL3-202 ENSG000001426… SH3BGRL3 Cont D N1
## # ℹ 17 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000270792.10 perc… D G 3 3 -0.222 2.45 0.842 0.842
## 2 ENST00000270792.10 perc… D I 3 3 2.16 2.22 0.15 0.3
## 3 ENST00000319041.6 perc… D G 3 3 -1.59 2.55 0.225 0.225
## 4 ENST00000319041.6 perc… D I 3 3 -3.60 2.08 0.065 0.131
## 5 ESPRESSO:chr1:910… perc… D G 3 3 1.50 2.27 0.259 0.259
## 6 ESPRESSO:chr1:910… perc… D I 3 3 8.74 2.25 0.009 0.018
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('RPS27')
## # A tibble: 27 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 Cont D N1
## 2 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 Cont D N2
## 3 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 Cont D N3
## 4 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… G N1
## 5 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… G N2
## 6 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… G N3
## 7 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… I N1
## 8 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… I N2
## 9 ENST00000477151.2 RPS27-202 ENSG0000017795… RPS27 siME… I N3
## 10 ENST00000493224.5 RPS27-203 ENSG0000017795… RPS27 Cont D N1
## # ℹ 17 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000651669.1 perce… D G 3 3 2.12 3.94 0.102 0.132
## 2 ENST00000651669.1 perce… D I 3 3 2.52 3.97 0.066 0.132
## 3 ENST00000477151.2 perce… D G 3 3 1.08 3.80 0.343 0.686
## 4 ENST00000477151.2 perce… D I 3 3 -0.180 3.81 0.866 0.866
## 5 ENST00000493224.5 perce… D G 3 3 -2.51 3.57 0.073 0.073
## 6 ENST00000493224.5 perce… D I 3 3 -3.71 3.23 0.03 0.06
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('ARPC1B')
## # A tibble: 45 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B Cont D N1
## 2 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B Cont D N2
## 3 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B Cont D N3
## 4 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… G N1
## 5 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… G N2
## 6 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… G N3
## 7 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… I N1
## 8 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… I N2
## 9 ENST00000458033.6 ARPC1B-212 ENSG0000013042… ARPC1B siME… I N3
## 10 ENST00000468337.2 ARPC1B-214 ENSG0000013042… ARPC1B Cont D N1
## # ℹ 35 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000646101.2 perc… D G 3 3 -1.64 3.77 0.181 0.312
## 2 ENST00000646101.2 perc… D I 3 3 -1.86 3.15 0.156 0.312
## 3 ENST00000481997.5 perc… D G 3 3 0.488 2.87 0.66 0.66
## 4 ENST00000481997.5 perc… D I 3 3 -1.44 2.52 0.261 0.522
## 5 ENST00000468337.2 perc… D G 3 3 -0.0789 2.27 0.944 0.944
## 6 ENST00000468337.2 perc… D I 3 3 -1.25 3.24 0.293 0.586
## 7 ENST00000645391.1 perc… D G 3 3 1.45 2.77 0.25 0.25
## 8 ENST00000645391.1 perc… D I 3 3 8.90 2.80 0.004 0.008
## 9 ENST00000458033.6 perc… D G 3 3 2.96 2.33 0.081 0.081
## 10 ENST00000458033.6 perc… D I 3 3 5.55 3.66 0.007 0.013
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('RPS3A')
## # A tibble: 18 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A Cont D N1
## 2 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A Cont D N2
## 3 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A Cont D N3
## 4 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… G N1
## 5 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… G N2
## 6 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… G N3
## 7 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… I N1
## 8 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… I N2
## 9 ENST00000274065.9 RPS3A-201 ENSG0000014542… RPS3A siME… I N3
## 10 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A Cont D N1
## 11 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A Cont D N2
## 12 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A Cont D N3
## 13 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… G N1
## 14 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… G N2
## 15 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… G N3
## 16 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… I N1
## 17 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… I N2
## 18 ENST00000507327.5 RPS3A-205 ENSG0000014542… RPS3A siME… I N3
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000274065.9 perce… D G 3 3 2.75 2.59 0.083 0.167
## 2 ENST00000274065.9 perce… D I 3 3 1.38 3.43 0.249 0.249
## 3 ENST00000507327.5 perce… D G 3 3 -0.699 3.74 0.526 0.526
## 4 ENST00000507327.5 perce… D I 3 3 1.93 3.38 0.138 0.276
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_percentage_isoforms('RPL41')
## # A tibble: 45 × 9
## transcript_id transcript_name gene_id gene_name type si rep
## <fct> <chr> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 Cont D N1
## 2 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 Cont D N2
## 3 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 Cont D N3
## 4 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… G N1
## 5 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… G N2
## 6 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… G N3
## 7 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… I N1
## 8 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… I N2
## 9 ENST00000501597.3 RPL41-202 ENSG0000022911… RPL41 siME… I N3
## 10 ENST00000546591.6 RPL41-204 ENSG0000022911… RPL41 Cont D N1
## # ℹ 35 more rows
## # ℹ 2 more variables: sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000546591.6 perc… D G 3 3 0.633 3.22 0.569 0.569
## 2 ENST00000546591.6 perc… D I 3 3 1.67 3.40 0.182 0.364
## 3 ENST00000501597.3 perc… D G 3 3 -2.06 3.87 0.11 0.11
## 4 ENST00000501597.3 perc… D I 3 3 -3.97 3.70 0.019 0.038
## 5 ENST00000546654.1 perc… D G 3 3 -0.576 2.04 0.622 0.74
## 6 ENST00000546654.1 perc… D I 3 3 1.14 2.03 0.37 0.74
## 7 ESPRESSO:chr12:4… perc… D G 3 3 1.56 2.71 0.225 0.45
## 8 ESPRESSO:chr12:4… perc… D I 3 3 -0.807 4.00 0.465 0.465
## 9 ENST00000552314.1 perc… D G 3 3 0.0147 2.22 0.989 0.989
## 10 ENST00000552314.1 perc… D I 3 3 -1.33 2.26 0.303 0.606
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('ERLIN2')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 36 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 2 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 3 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 4 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 5 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 6 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 7 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 8 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 9 ENST00000519638.3 protein_coding ERLIN2-204 ENSG00… protein_… ERLIN2
## 10 ENST00000648919.1 protein_coding ERLIN2-211 ENSG00… protein_… ERLIN2
## # ℹ 26 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 8 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000519638.3 norm… D G 3 3 -3.58 2.05 0.067 0.086
## 2 ENST00000519638.3 norm… D I 3 3 -4.54 2.05 0.043 0.086
## 3 ENST00000648919.1 norm… D G 3 3 5.22 3.48 0.009 0.019
## 4 ENST00000648919.1 norm… D I 3 3 1.88 3.95 0.135 0.135
## 5 ENST00000335171.10 norm… D G 3 3 -0.0642 2.04 0.955 0.955
## 6 ENST00000335171.10 norm… D I 3 3 16.8 2 0.004 0.007
## 7 ENST00000523887.5 norm… D G 3 3 1.67 2 0.237 0.474
## 8 ENST00000523887.5 norm… D I 3 3 0.287 3.53 0.79 0.79
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('NAP1L1')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 45 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 2 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 3 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 4 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 5 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 6 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 7 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 8 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 9 ENST00000549596.5 protein_coding NAP1L1-215 ENSG00… protein_… NAP1L1
## 10 ENST00000393263.7 protein_coding NAP1L1-201 ENSG00… protein_… NAP1L1
## # ℹ 35 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000549596.5 norm… D G 3 3 -3.08 3.37 0.046 0.08
## 2 ENST00000549596.5 norm… D I 3 3 -3.84 2.60 0.04 0.08
## 3 ENST00000393263.7 norm… D G 3 3 -0.00906 3.49 0.993 0.993
## 4 ENST00000393263.7 norm… D I 3 3 -3.44 2.55 0.053 0.105
## 5 ESPRESSO:chr12:4… norm… D G 3 3 -0.270 3.93 0.801 0.801
## 6 ESPRESSO:chr12:4… norm… D I 3 3 -1.10 2.32 0.371 0.742
## 7 ENST00000618691.5 norm… D G 3 3 8.93 2.01 0.012 0.012
## 8 ENST00000618691.5 norm… D I 3 3 15.5 2.02 0.004 0.008
## 9 ENST00000535020.6 norm… D G 3 3 -3.16 3.75 0.037 0.075
## 10 ENST00000535020.6 norm… D I 3 3 -0.438 2.47 0.697 0.697
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('RNF216')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 18 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 2 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 3 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 4 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 5 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 6 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 7 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 8 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 9 ENST00000425013.6 protein_coding RNF216-205 ENSG00… protein_… RNF216
## 10 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 11 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 12 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 13 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 14 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 15 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 16 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 17 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## 18 ENST00000389902.8 protein_coding RNF216-202 ENSG00… protein_… RNF216
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000425013.6 normc… D G 3 3 5.66 3.94 0.005 0.006
## 2 ENST00000425013.6 normc… D I 3 3 6.72 3.90 0.003 0.006
## 3 ENST00000389902.8 normc… D G 3 3 -1.79 2 0.215 0.43
## 4 ENST00000389902.8 normc… D I 3 3 -1.68 2 0.235 0.43
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('ITGB1-DT')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 18 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 2 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 3 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 4 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 5 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 6 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 7 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 8 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 9 ESPRESSO:chr10:2… <NA> <NA> ENSG00… lncRNA ITGB1-DT
## 10 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 11 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 12 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 13 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 14 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 15 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 16 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 17 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## 18 ENST00000691671.1 lncRNA ITGB1-DT-206 ENSG00… lncRNA ITGB1-DT
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ESPRESSO:chr10:24… norm… D G 3 3 -6.15 2 0.025 0.051
## 2 ESPRESSO:chr10:24… norm… D I 3 3 -1.46 2 0.282 0.282
## 3 ENST00000691671.1 norm… D G 3 3 -0.0903 3.99 0.932 0.932
## 4 ENST00000691671.1 norm… D I 3 3 1.02 3.63 0.372 0.744
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('DEPDC1')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 27 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 2 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 3 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 4 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 5 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 6 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 7 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 8 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 9 ENST00000370966.9 protein_coding DEPDC1-201 ENSG00… protein_… DEPDC1
## 10 ENST00000456315.7 protein_coding DEPDC1-202 ENSG00… protein_… DEPDC1
## # ℹ 17 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000370966.9 normc… D G 3 3 -3.47 2.36 0.058 0.116
## 2 ENST00000370966.9 normc… D I 3 3 -2.31 3.70 0.087 0.116
## 3 ENST00000456315.7 normc… D G 3 3 2.18 2.24 0.147 0.222
## 4 ENST00000456315.7 normc… D I 3 3 2.65 2.12 0.111 0.222
## 5 ENST00000488146.5 normc… D G 3 3 -1 2 0.423 0.846
## 6 ENST00000488146.5 normc… D I 3 3 -1 2 0.423 0.846
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('EZR')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 18 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 2 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 3 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 4 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 5 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 6 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 7 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 8 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 9 ENST00000367075.4 protein_coding EZR-202 ENSG00… protein_… EZR
## 10 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 11 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 12 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 13 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 14 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 15 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 16 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 17 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## 18 ENST00000337147.… protein_coding EZR-201 ENSG00… protein_… EZR
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000367075.4 norm… D G 3 3 1.76 3.94 0.155 0.268
## 2 ENST00000367075.4 norm… D I 3 3 1.88 3.93 0.134 0.268
## 3 ENST00000337147.11 norm… D G 3 3 -2.77 3.97 0.051 0.102
## 4 ENST00000337147.11 norm… D I 3 3 -0.392 3.43 0.718 0.718
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('ITGA6')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 27 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 2 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 3 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 4 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 5 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 6 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 7 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 8 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 9 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 10 ENST00000409080.6 protein_coding ITGA6-202 ENSG00… protein_… ITGA6
## # ℹ 17 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000684293.1 norm… D G 3 3 -11.5 2.69 0.002 0.005
## 2 ENST00000684293.1 norm… D I 3 3 -6.43 3.14 0.007 0.007
## 3 ENST00000409080.6 norm… D G 3 3 2.12 2.37 0.147 0.294
## 4 ENST00000409080.6 norm… D I 3 3 1.71 2.18 0.219 0.294
## 5 ENST00000264107.12 norm… D G 3 3 -3.87 3.68 0.021 0.042
## 6 ENST00000264107.12 norm… D I 3 3 -1.11 2.62 0.357 0.357
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('UAP1')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 36 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 2 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 3 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 4 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 5 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 6 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 7 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 8 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 9 ENST00000367926.9 protein_coding UAP1-202 ENSG00… protein_… UAP1
## 10 ENST00000367925.6 protein_coding UAP1-201 ENSG00… protein_… UAP1
## # ℹ 26 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 8 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000367926.9 norm… D G 3 3 4.65 2.04 0.041 0.083
## 2 ENST00000367926.9 norm… D I 3 3 0.880 3.57 0.434 0.434
## 3 ENST00000367925.6 norm… D G 3 3 -1.63 3.35 0.192 0.384
## 4 ENST00000367925.6 norm… D I 3 3 -0.607 3.78 0.578 0.578
## 5 ESPRESSO:chr1:173… norm… D G 3 3 -0.362 3.41 0.738 0.986
## 6 ESPRESSO:chr1:173… norm… D I 3 3 0.757 3.86 0.493 0.986
## 7 ENST00000486089.1 norm… D G 3 3 -0.689 4.00 0.528 1
## 8 ENST00000486089.1 norm… D I 3 3 -0.595 3.99 0.584 1
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('ITGA6')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 27 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 2 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 3 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 4 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 5 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 6 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 7 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 8 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 9 ENST00000684293.1 protein_coding ITGA6-212 ENSG00… protein_… ITGA6
## 10 ENST00000409080.6 protein_coding ITGA6-202 ENSG00… protein_… ITGA6
## # ℹ 17 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 6 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000684293.1 norm… D G 3 3 -11.5 2.69 0.002 0.005
## 2 ENST00000684293.1 norm… D I 3 3 -6.43 3.14 0.007 0.007
## 3 ENST00000409080.6 norm… D G 3 3 2.12 2.37 0.147 0.294
## 4 ENST00000409080.6 norm… D I 3 3 1.71 2.18 0.219 0.294
## 5 ENST00000264107.12 norm… D G 3 3 -3.87 3.68 0.021 0.042
## 6 ENST00000264107.12 norm… D I 3 3 -1.11 2.62 0.357 0.357
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('ADD3')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 36 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 2 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 3 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 4 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 5 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 6 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 7 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 8 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 9 ENST00000356080.9 protein_coding ADD3-202 ENSG00… protein_… ADD3
## 10 ENST00000277900.… protein_coding ADD3-201 ENSG00… protein_… ADD3
## # ℹ 26 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 8 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000356080.9 norm… D G 3 3 7.84 4.00 0.001 0.003
## 2 ENST00000356080.9 norm… D I 3 3 1.33 2.86 0.281 0.281
## 3 ENST00000277900.12 norm… D G 3 3 -3.69 3.89 0.022 0.027
## 4 ENST00000277900.12 norm… D I 3 3 -6.27 2.51 0.014 0.027
## 5 ENST00000488799.5 norm… D G 3 3 -0.877 3.87 0.432 0.432
## 6 ENST00000488799.5 norm… D I 3 3 -1.77 2.10 0.212 0.424
## 7 ENST00000360162.7 norm… D G 3 3 -2.39 2.19 0.128 0.128
## 8 ENST00000360162.7 norm… D I 3 3 -2.66 3.64 0.062 0.124
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('SAP18')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 45 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 2 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 3 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 4 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 5 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 6 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 7 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 8 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 9 ENST00000607003.5 protein_coding SAP18-207 ENSG00… protein_… SAP18
## 10 ENST00000382533.9 protein_coding SAP18-201 ENSG00… protein_… SAP18
## # ℹ 35 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000607… norm… D G 3 3 -14.8 3.23 4.42e-4 8.84e-4
## 2 ENST00000607… norm… D I 3 3 1.27 3.31 2.86e-1 2.86e-1
## 3 ENST00000382… norm… D G 3 3 16.5 2 4 e-3 7 e-3
## 4 ENST00000382… norm… D I 3 3 1.27 3.31 2.86e-1 2.86e-1
## 5 ENST00000621… norm… D G 3 3 16.5 2 4 e-3 7 e-3
## 6 ENST00000621… norm… D I 3 3 1.27 3.31 2.86e-1 2.86e-1
## 7 ENST00000450… norm… D G 3 3 1.37 2.23 2.93e-1 5.86e-1
## 8 ENST00000450… norm… D I 3 3 0.847 3.40 4.53e-1 5.86e-1
## 9 ENST00000471… norm… D G 3 3 1.95 2.76 1.55e-1 3.1 e-1
## 10 ENST00000471… norm… D I 3 3 1.25 3.16 2.97e-1 3.1 e-1
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('TSPAN8')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 36 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 2 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 3 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 4 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 5 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 6 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 7 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 8 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 9 ENST00000546561.2 protein_coding TSPAN8-203 ENSG00… protein_… TSPAN8
## 10 ENST00000247829.8 protein_coding TSPAN8-201 ENSG00… protein_… TSPAN8
## # ℹ 26 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 8 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000546561.2 norm… D G 3 3 6.21 3.95 0.004 0.007
## 2 ENST00000546561.2 norm… D I 3 3 7.43 3.05 0.005 0.007
## 3 ENST00000247829.8 norm… D G 3 3 -5.07 3.30 0.012 0.023
## 4 ENST00000247829.8 norm… D I 3 3 -2.83 2.12 0.099 0.099
## 5 ESPRESSO:chr12:41… norm… D G 3 3 -1.62 2.56 0.219 0.219
## 6 ESPRESSO:chr12:41… norm… D I 3 3 -2.76 3.23 0.065 0.13
## 7 ENST00000552128.2 norm… D G 3 3 0.220 2.73 0.842 0.842
## 8 ENST00000552128.2 norm… D I 3 3 2.09 3.80 0.109 0.218
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('PFDN5')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 54 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 2 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 3 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 4 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 5 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 6 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 7 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 8 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 9 ENST00000334478.9 protein_coding PFDN5-202 ENSG00… protein_… PFDN5
## 10 ENST00000351500.7 protein_coding PFDN5-203 ENSG00… protein_… PFDN5
## # ℹ 44 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 12 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000334478.9 norm… D G 3 3 -6.55 3.98 0.003 0.006
## 2 ENST00000334478.9 norm… D I 3 3 -1.08 2.10 0.388 0.388
## 3 ENST00000351500.7 norm… D G 3 3 -7.40 2.08 0.016 0.025
## 4 ENST00000351500.7 norm… D I 3 3 -8.34 2.07 0.013 0.025
## 5 ENST00000243040.… norm… D G 3 3 0.374 2.44 0.738 0.738
## 6 ENST00000243040.… norm… D I 3 3 -3.63 3.25 0.031 0.063
## 7 ENST00000548984.1 norm… D G 3 3 -1.02 3.99 0.365 0.73
## 8 ENST00000548984.1 norm… D I 3 3 -1.09 2.54 0.368 0.73
## 9 ENST00000550846.5 norm… D G 3 3 -2.17 2.63 0.131 0.131
## 10 ENST00000550846.5 norm… D I 3 3 -4.05 3.35 0.022 0.044
## 11 ENST00000547130.6 norm… D G 3 3 -2.00 2.38 0.163 0.163
## 12 ENST00000547130.6 norm… D I 3 3 -7.15 3.13 0.005 0.01
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('RTN4')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 45 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 2 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 3 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 4 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 5 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 6 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 7 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 8 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 9 ENST00000394609.6 protein_coding RTN4-205 ENSG00… protein_… RTN4
## 10 ENST00000317610.… protein_coding RTN4-201 ENSG00… protein_… RTN4
## # ℹ 35 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 10 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000394… norm… D G 3 3 -0.674 3.15 5.46e-1 5.46e-1
## 2 ENST00000394… norm… D I 3 3 -20.1 3.82 5.11e-5 1.02e-4
## 3 ENST00000317… norm… D G 3 3 -3.22 3.49 3.9 e-2 7.8 e-2
## 4 ENST00000317… norm… D I 3 3 2.57 3.95 6.3 e-2 7.8 e-2
## 5 ENST00000491… norm… D G 3 3 0.411 3.24 7.07e-1 7.07e-1
## 6 ENST00000491… norm… D I 3 3 -1.75 3.98 1.55e-1 3.1 e-1
## 7 ENST00000485… norm… D G 3 3 -0.907 2.05 4.58e-1 5.5 e-1
## 8 ENST00000485… norm… D I 3 3 -1.49 2.00 2.75e-1 5.5 e-1
## 9 ENST00000357… norm… D G 3 3 -2.53 3.99 6.5 e-2 1.3 e-1
## 10 ENST00000357… norm… D I 3 3 -1.05 3.33 3.63e-1 3.63e-1
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('S100A4')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 54 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 2 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 3 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 4 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 5 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 6 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 7 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 8 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 9 ENST00000368716.9 protein_coding S100A4-204 ENSG00… protein_… S100A4
## 10 ENST00000481009.1 protein_coding… S100A4-206 ENSG00… protein_… S100A4
## # ℹ 44 more rows
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 12 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST00000368… norm… D G 3 3 -15.7 2.48 2 e-3 2 e-3
## 2 ENST00000368… norm… D I 3 3 -16.6 2.89 5.97e-4 1 e-3
## 3 ENST00000481… norm… D G 3 3 -11.4 3.57 6.23e-4 1 e-3
## 4 ENST00000481… norm… D I 3 3 -9.86 3.36 1 e-3 1 e-3
## 5 ENST00000354… norm… D G 3 3 52.9 2.59 5.29e-5 5.29e-5
## 6 ENST00000354… norm… D I 3 3 36.9 3.47 1.26e-5 2.52e-5
## 7 ENST00000468… norm… D G 3 3 0.0207 2.69 9.85e-1 9.85e-1
## 8 ENST00000468… norm… D I 3 3 -35.0 4.00 3.99e-6 7.98e-6
## 9 ENST00000368… norm… D G 3 3 -10.3 3.54 9.29e-4 2 e-3
## 10 ENST00000368… norm… D I 3 3 -6.85 3.48 4 e-3 4 e-3
## 11 ESPRESSO:chr… norm… D G 3 3 7.75 2.03 1.6 e-2 1.6 e-2
## 12 ESPRESSO:chr… norm… D I 3 3 11.1 2.88 2 e-3 4 e-3
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
barplot_normcount_isoforms('RPS24')
## Joining with `by = join_by(transcript_id, transcript_name)`
## # A tibble: 18 × 13
## transcript_id transcript_type transcript_name gene_id gene_type gene_name
## <fct> <chr> <chr> <chr> <chr> <chr>
## 1 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 2 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 3 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 4 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 5 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 6 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 7 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 8 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 9 ENST00000613865.5 protein_coding RPS24-214 ENSG00… protein_… RPS24
## 10 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 11 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 12 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 13 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 14 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 15 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 16 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 17 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## 18 ENST00000372360.9 protein_coding RPS24-202 ENSG00… protein_… RPS24
## # ℹ 7 more variables: seqname <chr>, type <chr>, si <chr>, rep <chr>,
## # normcount <dbl>, sum_normcount <dbl>, percent_transcript <dbl>
## # A tibble: 4 × 13
## transcript_id .y. group1 group2 n1 n2 statistic df p p.adj
## <fct> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl>
## 1 ENST000006138… norm… D G 3 3 -0.763 3.89 4.89e-1 4.89e-1
## 2 ENST000006138… norm… D I 3 3 -3.05 2.52 6.9 e-2 1.38e-1
## 3 ENST000003723… norm… D G 3 3 -53.3 3.57 2.59e-6 5.18e-6
## 4 ENST000003723… norm… D I 3 3 -6.14 2.09 2.3 e-2 2.3 e-2
## # ℹ 3 more variables: p.adj.signif <chr>, y.position <dbl>, groups <named list>
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
## No summary function supplied, defaulting to `mean_se()`
